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ABSTRACT 

The magnetic, chemically peculiar A-star HR 710 (HD 15144, ADS 1849A, AB Cet) 
has a close companion in a nearly circular orbit. Its 3-day period is unusually short 
for such stars. The system emits moderately hard x-rays, which likely come from a 
white dwarf secondary (ADS 1849 Ab). The Srii spectrum is very strong, and the 
resonance lines show similar core-nib structure to the stronger Ca K line. We place 
only loose constraints on a model. There are indications of a lower electron/gas pres- 
sure than expected from the star's parallax and brightness. Strong-line profiles and 
anomalous excitation/ionization indicate significant deviations from traditional mod- 
els, fixed by T e //, log<?, and abundances. Weak spectral-line profiles and wavelength 
shifts probably indicate abundance patches rather than the presence of a secondary 
(ADS 1849Ab) spectrum. Our UVES spectra are from only one epoch. Additional low- 
noise, high-resolution spectra are needed. We discuss the spectrum within the context 
of abundance stratification. 
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1 INTRODUCTION 

\ Roman (1949) listed HR 710 (HD 15144) among the prob- 
able members of the Ursa Major Stream. She noted the 
strengths of AA4077 and 4215. More recently, King, et al. 
(2003) enter a question mark in the column for "Final 
Memb." for HD 15144A. Tokovinin (1997) concludes that 
the system ADS 1849AB is a physical double, rather than 
a common proper motion pair, but the period must be very 
long. The primary, ADS 1849A (or is itself a spectroscopic 
binary (Aa and Ab) with a well-determined period of 2.998 
days (Bonsack 1981, Tokovinin 1997). This a very short pe- 
riod for a magnetic Ap star. One system, HD 200405 (cf. 
Carrier, et al. 2002), with an even shorter period is known- 
1.64 days. The HR 710 system appears to be single lined, so 
there has been a large uncertainty in the mass and luminos- 
ity of the companion. 

Htinsch, Schmitt, and Voges (1998, henceforth HSV) 
list HR 710 in their ROSAT survey of main sequence and 
subgiant stars of the Bright Star Catalogue. They found an 
x-ray luminosity of 4.2 x 10 29 erg s . This is among the top 

* Based on observations obtained at the European Southern Ob- 
servatory, Paranal and La Silla, Chile (ESO programmes 076. D- 
0169(A) and 076.C-0172(A)). 

f E-mail: cowley@umich.edu (CC); shubrig@eso.org (SH) 



third of the group. The hardness ratio they report, 0.12, is 
in the upper 15%. Thus HR 710 is significant, though not 
outstanding among these x-ray emitters. 

Czesla and Schmitt (2007) report Chandra observations 
of HR 710 as well as the companion ADS 1849B in the energy 
band 0.5 to 2keV Their x-ray flux for the primary is in 
reasonable agreement with the .ROSAT data. 

Extensive magnetic field measurements were made by 
Babcock (1958), and Bonsack (1981). The field varied be- 
tween —0.45 and —1.10 kG — it did not reverse. Bonsack's 
plot of H B s folded on the 3-day orbital period appears ran- 
dom. These results are consistent with nearly magnetic and 
rotational pole-on observations. Bonsack sought other peri- 
ods, and found marginal evidence for a 15.88-day period. He 
concluded < H c g >~ — 700 G, with an rms error of a single 
observation of 140 G. 

One of us, CRC, measured wavelengths of a 2.4 A mm -1 
spectrum from the Dominion Astrophysical Observatory 
in 1979. Contemporary notes include a remarks on the 
strengths of Srii and Can K. Additionally, they note the 
subordinate Srii line, 4305.45 "has huge wings," (see Fig. [1). 
Our description is similar to that of Babcock (1958), who 
also pointed out that Lines of Srii, like the K line, have 
strong, sharp cores with wide shallow wings. Much of the 
present study follows this lead. 
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The upper spectrum is of HR 710, while the lower is of 7 Equ. 
The latter spectrum is displaced downward by 0.3 for display 
purposes. Although both stars are classified as strontium stars, 
the wings on A4305 are arguably more pronounced in HR 710 
than in the sharper-lined spectrum of 7 Equ. 

Otherwise, the spectrum was not found unusual among 
the variegated spectra of CP types. Iron-group spectra, in- 
cluding Sell (but not Vn), are rich, while Yn and Zrll are 
not as strong as one might expect given the outstanding 
strength of Srn. The lanthanides are present, but not re- 
markable richness or in strength. The strong Eull lines at 
AA4129 and 4205, for example, are exceeded in strength by 
nearby lines of iron-peak elements. 

We know of no modern abundance study of HR 710. 
It was among several Ap stars analyzed in the 1960's, with 
rather crude methods by today's standards, and only se- 
lected elements were studied. For example, the work of 
Searle, Lungershausen, and Sargent (1966) was restricted 
to Ti, Mn, Cr, and Fe. We would agree with these authors 
that these particular elements are all enhanced by similar 
factors. 

HR 710 was among the stars studied by Cowley, et 
al. (2007, henceforth, CHCGW) for the calcium isotopic 
anomaly discovered by Castelli and Hubrig (2004). In the 
course of this work we confirmed that conventional models 
failed to account for the profiles of the strong Ca 11 lines. This 
was in complete agreement with results of Ryabchikova and 
her coworkers. See, for example, Ryabchikova, et al. (2002, 
Fig. 3). These workers have used a stratified-abundance 
model to fit Ca I and 11 and other metallic lines, as described 
in subsequent studies on numerous magnetic Ap stars (cf. 
Ryabchikova, Kochukhov, and Bagnulo 2007). 

For 7 Equ, we found that essentially the same model 
fit the A8542 line of the infrared triplet (henceforth, IRT, 
cf. Fig. 7 of CHCGW). The same stratification provided a 
reasonably good fit to the A8498-line. 

In the course of research leading to CHCGW a number 
of lines were synthesized that were not discussed in that 
work. Provisional calculations for HR 710 showed we would 
have difficulty matching the profiles of the IRT with the 
same stratification that fit the Call K-line. 

The unusual strength of Sr 11 in HR 710 deserves closer 
study. Would a stratified model be required to match the 
profiles of the Srn resonance lines AA4077 and 4215, or the 



strong subordinate lines? If so, would the study of these fea- 
tures add to our understanding of the structure of the pho- 
tospheres of magnetic stars? The sections below deal with 
these questions. 

We can place only loose constraints on the physical con- 
ditions in the atmosphere of HR 710. In addition to line 
core/wing and ionization and excitation anomalies, there is 
evidence for abundance patches on the surface, which can 
be seen in some of the line profiles. 



2 OBSERVATIONS 

The ESO spectra used in this study were obtained on 19 
September 2005 with the UV- Visual Echelle Spectrograph 
(UVES) at UT2 of the VLT. The spectra used the new 
(November 2004) standard Dichroic settings were described 
by CHCGW, covering 3290 to 9460 A. This configuration 
makes it possible to measure the A8542-line, which had pre- 
viously been unavailable because of an order gap. The re- 
solving power is 80 000 in the blue-violet, and 110 000 in the 
red. The S/N ratios of the spectra differ, but range from 
about 100 to 400. 



3 REDUCTIONS 

The UVES spectra were reduced by the UVES pipeline Data 
Reduction Software (version 2.5), which is an evolved version 
of the echelle context of MIDAS. The spectra were mildly 
Fourier filtered before the wavelengths were measured at 
Michigan. 

We choose the continuum by picking high points of the 
relative flux using console displays 40 A wide. In some cases, 
a Paschen line or water- vapor band would take up most of 
the display, and no continuum point would be chosen. The 
high points were connected by a spline and by linear interpo- 
lation. After some experimentation the linear fit was given 
weight 2, the spline 1, and the weighted average adopted. 
The continuum height is surely uncertain to a few per cent. 
Nevertheless we note that in the illustrations that follow, no 
additional (ad hoc) adjustments were made to the continua. 

Wavelength measurements from 4527 to 4652 A are 
from 2.4 A mm -1 spectra (12766U, and V) obtained at the 
Dominion Astrophysical Observatory in 1979. The instru- 
ment has been described by Fletcher et al. (1980). 



4 THE ATMOSPHERE OF HR 710 
4.1 Preliminary remarks 

Ap stars were once thought to have abnormal abundances, 
but normal atmospheres. We now know that this is not the 
case. The profiles of the low Balmer and the strong metal- 
lic lines cannot be reproduced by standard plane-parallel 
atmospheres. Simple adjustments to temperature, gravity, 
and composition do not allow one to resolve inconsistencies 
in ionization and excitation. These difficulties have become 
apparent relatively recently. Horizontal chemical inhomo- 
geneities (abundance patches) have been known in magnetic 
stars for many decades. 

Extensive observations by Babcock (1958) and Bonsack 
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(1981) show that spectral variations in HR 710 are mild or 
absent. Nevertheless, abundance patches apparently influ- 
ence the spectrum. We acknowledge this, among other com- 
plications, but do not treat it in the present study. 

4.2 Photometry 

We used several codes based on Stromgren or Geneva pho- 
tometry to try to fix T e // an d logg. In addition, we consid- 
ered ionization balance, Balmer, and Paschen profiles. The 
various implementations give a fairly wide range of values 
for Teff (8400 to 8600 K) and logg (3.0 to 4.0). 

An implementation of Geneva Photometry by Kaiser 

(2006) , based on Kiinzli, et al. (1997) includes an adjust- 
ment for [Fe/H]. The code gives values of logg that de- 
crease with increasing [Fe/H]. Allowed entries for [Fe/H] are 
— 1, 0, and +1. The relation for both logg and T e // is very 
nearly linear. We used photometric measurements for HR 
710 from the Simbad data base to obtain second-order fits, 
with very small quadratic terms. With [Fe/H] = +1.0, we 
obtain T e // = 8460 K and logg = 3.6. This somewhat low 
gravity is supported by our ionization balance of neutral 
and singly ionized Fe, where we get best agreement with 
logg = 3.0. This ionization balance (see below) gives a tem- 
perature nearly 200 K higher. 

Codes based on Stromgren photometry give signifi- 
cantly higher gravities, though they do not, to our knowl- 
edge, include a provision to allow for metallicity. 

The higher surface gravities agree with Hubrig, et al. 

(2007) , who give logg = 4.2. The temperature estimate from 
this study is 8433 K. 

4.3 Balmer lines 

Hydrogen line strengths depend sensitively on temperature 
in cooler stars, and on surface gravity in hotter ones. HR 710 
is at the high-temperature end of the region where the line 
strengths change from primarily temperature to primarily 
gravity sensitivity. 

All Balmer calculations are based on the profiles from 
Stehle & Hutcheon (1999). 

Fig. [2] shows the observed and calculated profile of Ha 
for models with T eff = 8600 K, and logy = 3.0 and 4.0. The 
fit is arguably better with the lower-gravity model (upper 
smooth curve). Metallic lines have been omitted from the 
calculation to avoid confusion. The tendency of the lower 
Balmer cores of cooler magnetic Ap stars to be deeper, 
and sharper than the observations was called the core-wing 
anomaly (Cowley, et al. 2001). LTE calculations of the Ha 
core are never as deep as the observations. In the case of 
HR 710, departures of the calculations from the observed 
profiles in the regions outside the inner core are of interest. 
For Ha and especially for H7 and H<5 the calculations do not 
fit the observations until the far wings are reached. 

The case for the lower gravity is somewhat stronger with 
H7 (Fig. E| and H<5 (not shown). 

Fig. [4] shows the region of the Balmer confluence. H18 
(3691.56) is surely present, and H19 (3686.83) possibly so, 
though the feature is clearly strengthened by blending pri- 
marily with Cm 3687.03 and Fe 1-75 3687.10. The Multiplet 
Number (75) is from Moore(1945). The strong Cr n line does 
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Figure 2. Ha profile in HR 710. The calculated profiles from 
models with logg = 3.0 (upper curve), and 4.0 (lower curve). In 
both cases, the assumed temperature is 8600 K. 
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Figure 3. H7 profile in HR 710. The calculated profiles from 
models with logg = 3.0 (upper curve), and 4.0 (lower curve). In 
both cases, the assumed temperature is 8600 K. 



not appear in the Multiplet Tables. The Inglis- Teller rela- 
tion (cf. Cowley 1970, Eq. AI-5) yields log N e = 13.4 if the 
last Balmer line is H19. This electron density is reached in 
our models with T e ff = 8600 K, and with logg = 4 at a 
logarithmic optical depth logT3g4o = —0.7. For the logg = 3 
model the corresponding optical depth is —0.08. 

We have used the optical depth at A3640, short of the 
Balmer Jump. This should closely approximate the opacity 
where the confluence occurs, including the opacity due to 
overlapping higher Balmer lines. The basis of this statement 
may be found, among other places, in Cowley (2000, see also 
Appendix A). 

We expect the electron density obtained by the Inglis- 
Teller formula would correspond to N e at optical depth unity 
for the wavelength of the Balmer jump. Thus, the last count- 
able Balmer line appears to give mild support for a lower- 
log g model. However a plot similar to Fig.|4]for 7 Equ shows 
Balmer components possibly as high as H20. Moreover, in 
two other main-sequence stars, v Cap (B9.5 V) and a Boo 
(F2 V), we can count to H22 or H23 and H19 or 20 respec- 
tively. A detailed study of the usefulness of the last count- 
able Balmer line as an indicator of surface gravity would be 
worth while. 
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Figure 5. The Paschen line P14. The smooth, solid curve was 
calculated using a model with T e ff = 8600, and logg = 4. 
The dashed line shows the profile from the lower-gravity model, 
logfl = 3. 



Figure 4. Region of the Balmer confluence. The spectrum has 
been rather severely filtered in an attempt to make the Balmer 
lines stand out. 



4.4 Paschen lines 

The lower Paschen members are beyond our wavelength cov- 
erage. The first usable line is A9229, P9. The lowest Paschen 
members, P9, P10, and Pll, are clearly fit better with a 
log g = 4 model than log g = 3. All Paschen line calculations 
used Lemke's (1997) revised VCS profiles (Vidal, Cooper & 
Smith 1973), because the Stehle-Hutcheon tables do not ex- 
tend to large enough values of N e for the higher Paschen 
lines. 

At A8750, P12, the lower-gravity model provides a bet- 
ter fit to the wings than the log g = 4 model, but the core is 
too deep. If we were to lower the observed continuum by 2 to 
3%, the overall fit, core and wings is significantly better with 
the higher-gravity model. Fig. [5] shows the Paschen line at 
A8598, P14, which falls between the longer-wavelength pair 
of the IRT. Again, the higher- gravity model provides a bet- 
ter fit. 

One may surely count Paschen lines to P18. It is ar- 
guable whether higher Paschen members are present. Thus 
the limiting Balmer and Paschen series are in agreement 
with one another, but a definitive interpretation awaits fur- 
ther study. Better observational material, specifically taken 
and reduced for the purpose of investigation of series con- 
fluence would also be useful. 

The Paschen lines P13, P14, P15, and P16 are blended 
with present special problems because they are near the se- 
ries confluence. Details are explained in Appendix A. 



5 IONIZATION EQUILIBRIUM 

Lines of Fel and II were chosen from a compilation of os- 
cillator strengths by Fuhr and Wiese (2006). We used our 
line identification list (see URL below) to exclude cases with 
obvious blends. We then calculated abundances for models 
with logg = 3 and 4 for T eff ranging from 8000 to 8800 K. 
Calculations for all of these cases were made with, and with- 



out stratification. The stratification assumed was based on 
an analytic formula: 



b) 



a ■ \x + d\ 2 



(1) 



In any given layer, the abundance of an element is mul- 
tiplied by this g(x). Here, x — log(r5ooo), and the negative 
sign is taken for x ^ — d. This function is a smoothed step 
of height 1 — b, with the center of the jump located at depth 
logT5ooo = —d. At large depths, g(x) is nearly unity, while 
for the smallest depths, it has the value b. The Gaussian 1/e 
width of the step is 1/ yfa. 

For all spectra other than those of Call and Srll, "stan- 
dard parameters" were used: a = 6.7, log b = —4, and d = 1. 
We cannot claim any special validity for these particular val- 
ues; they were used in CHCGW to reconcile profiles of the 
Ca K-line and one line of the IRT. We have already explored 
a wide domain of parameter space, and must leave possible 
adjustments of the stratification parameters to account for 
anomalous ionization and excitation to future work. 

Tables of wavelengths and equivalent widths may be 
found on our web page: 

http : //www. astro . lsa.umich. edu/~cowley/hr710 

For iron, abundances were found for 5 temperatures and 
2 gravities, above and below the Balmer Jump, and with 
and without stratification using the standard parameters. 
We made plots of abundances versus equivalent widths and 
excitation potential, and used the former plots to fix the 
microturbulent parameter. Trends for iron were removed by 
the value £t = 3.5. This value is somewhat large, but it also 
must account for the Zeeman broadening, which was not 
large enough to yield resolved patterns. 

The results are illustrated in Fig. [6] Here, and elsewhere 
in this paper, 'N to t' means the sum of relative abundances 
of all elements. If log(H) = 12, then N to t = 1-03 • 10 12 . We 
used 70 lines of Fe I and 32 lines of Fe II, all with wavelengths 
longer than 3650 A. This plot is the primary basis for con- 
sideration of low-gravity (or low pressure) models. 

For the 8600 K, logg = 3 model, with no stratification, 
we find 



log (Fe/Nt. 



-3.93 (Fel); log (Fe/N t 



-3.89 (Fell) 
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Figure 6. Ionization equilibrium: Fel and II. The dashed lines 
are Fe I (upper) and Fe II (lower) using the standard stratification 
parameters. The lower lines are Fe I (mostly upper) and Fell 
(mostly lower) without stratification. Values of log(Fe/Ntot) refer 
to the deep photosphere (prior to the application of the function 
g(x)). The four curves are based on lines with wavelengths longer 
than 3650 A. 



For iron, we made separate abundance determinations 
for lines below the Balmer Jump. For 43 Fel lines, we ob- 
tained log(Fe/Ntot) — —3.86. This is less than 0.1 dex differ- 
ent from the value found for lines above the Balmer Jump. 
Only 8 lines of Fell were found usable below the Balmer 
Jump. These yielded log(Fe/N to t) = —4.02. In view of the 
smaller number of lines, we do not consider this value sig- 
nificantly different from the others. All calculations of iron 
below the Balmer Jump were made without assuming strat- 
ification. We see no compelling evidence for stratification in 
iron. 

The next richest spectra to Fe I and n belong to Cr I and 
II. We analyzed 36 Crl lines and 47 Cm lines. Oscillator 
strengths were taken from VALD (see Kupka, et al. 1999), 
but we omitted lines from intersystem transitions, where the 
multiplicity of the upper and lower terms are different. Ex- 
perience has shown that ab initio and semi-empirical calcu- 
lations for LS-permitted transitions are somewhat more reli- 
able than those for intersystem lines. Reasonable agreement 
was found for abundances from the Crl and II lines with 
wavelengths longer than 3650 A, from the 8600 K, logg = 3 
model, with no stratification: 

log(Cr/N t ot) = -4.00 (Crl); log(Cr/N to t) = -4.14 (Cr II) 

Several other elements are observed in two stages of ion- 
ization, but our results from them were less concordant. In 
particular, the spectra of Til and n may be brought into 
better agreement by adjusting model parameters or intro- 
ducing stratification. Details are tedious, and are omitted, 
pending a later investigation. 



6 CALCIUM 

6.1 The weaker Cal and II lines 

We measured equivalent widths of 23 lines of Ca I, all with 
wavelengths longer than the Balmer Jump. The equivalent 
widths spread evenly between 10 and 112 mA. With the 
8600 K, logg = 3 model, we obtained a mean log(Ca/Ntot) 
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Figure 7. The Can K-line, observed (solid gray) and calcu- 
lated with stratification. The stratification parameters for g(x) 
are: a = 6.7, logfe = -7.0, d = 0.25, and log (Ca/Ntot) = -3.10. 
The heavy dashed line is for log (Ca/Ntot) = —4.65, and no strat- 
ification. 



of —4.65 ±0.08 (£t = 3.5). There is no significant trend of 
the abundance with either equivalent width or excitation 
potential; the latter span a narrow range of 1.88 to 2.93 eV. 

The Ca I resonance line, A4227 was not included in the 
above average. It is only moderately strong (W\ — 169 mA), 
and its core is in not at all unusual, like those of the K- 
line or AA4077 or 4215 of SrII. We can fit the profile with 
log(Ca/N to t) of —4.65, and lowering the £ t to 1.5 kms -1 . We 
conclude A4227 is in reasonable agreement with the weaker 
Cal lines. 

We found only 7 useful Can lines, exclusive of H and K, 
and the IRT. The mean value of log(Ca/N to t) from all 7 lines 
was -4.99 ±0.33, with the 8600 K, logg = 3 model. With 
these lines there were systematic trends both with excitation 
potential and equivalent width. Four lines, with excitation 
potentials averaging 7.2 eV gave log(Ca/Ntot) = —5.5, while 
three lines with average Xcx of 9.0 eV gave log(Ca/Ntot) = 
—4.2. The tendency of high-excitation lines to give anoma- 
lously high abundance is a characteristic of the newly dis- 
cussed anomalies of CP stars (cf. Ryabchikova, et al. 2004). 
However, the abundance averages from Cal and II are within 
the standard errors; we have not attempted to deal with the 
systematic trends. 

One Call line is below the Balmer Jump: 3461.87 (x = 
7.51). The abundance from this single line is log(Ca/N to t) = 
—4.6. We make no judgement about stratification from this 
single line. 



6.2 The K-line 

The Can K-line exhibits the full symptoms of stratification. 
An attempt to fit it with the abundance obtained from the 
Cal or Can lines (cf. above) results in a deep, broad core, 
as shown in Fig. [7] The plot is quite similar to Fig. 3 of 
Ryabchikova, et al. (2002) for the star 7 Equ. Here, and in 
similar figures to follow, only token efforts were made to fit 
the atomic line spectrum apart from the feature of direct 
interest. 

The stratification parameters given in the caption to 
Fig. [7] were determined by trial and error. They imply an 
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Figure 8. The Call A8498-line, observed (solid gray) and cal- 
culated (solid black) with stratification. The stratification pa- 
rameters for g(x) are: a = 6.7, logb = —5.70, d = 0.25. 
log(Ca/Ntot) = —2.12. The dashed line is for the parameters used 
to fit the Can K-line (see Fig. 0. Both calculations employed a 
model with log g = 4. 
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Figure 9. The Can A8542-line, observed (solid gray) and calcu- 
lated (solid black) with stratification. The stratification parame- 
ters for g(x) are: a = 1.0, logb = —6.70, d = —1.0. log(Ca/Ntot = 
— 1.52. The function g(x) differs from unity throughout the at- 
mosphere, so the effective log(Ca/N to t) at log(T5ooo) = 0.0 is 
roughly —2.52, in fair agreement with the value necessary to fit 
the A8498 line. The dashed line is for the parameters used to fit 
the A8498 line. Both calculations use logg = 4. 



abundance jump of 7 orders of magnitude between optical 
depths log (r 5 ooo) = —0.7 and +0.3. 



6.3 The IRT 

Lines of the IRT manifest stratification in a way that is 
somewhat similar to the Can K-line, though the total ab- 
sorption is much less, and the central absorption is more sig- 
nificant relative to the total profile. We must keep in mind 
that the possible influence of surficial abundance patches 
which could have an influence on the profiles of the IRT. In 
view of the better fits with log g — 4 for P14, which falls be- 
tween A8498 and A8542 of the IRT, as well as lowerPaschen 
all calculations illustrated for the IRT have used the higher 
gravity. 

Fig [5] shows observed and computed profiles for A8498 
in HR 710. Damping wings are present, but the stellar pro- 
files change rapidly in structure, from core to wing. These 
shapes can arise if the abundance of calcium varies with a 
stratification function g(x) as described above. 

Note the extraordinarily high value of calcium required 
for the deep photosphere, and the large abundance gradient 
implied by the parameters of the function g(x). 

The fit to the strongest line of the IRT, A8542, is 
shown in Fig. [9] The wings of this line are significantly 
stronger than those of the weaker A8498 line. The best 
fit (solid line) used different stratification parameters from 
those used to fit the weaker A8498 line (dashed line), though 
the abundance in the deep photosphere is roughly similar: 
log(Ca/Ntot) = —2.5 vs. —2.1. Several lines of neutral cal- 
cium are in the region of Fig. [5] Two of the lines, at 8525.72 
and 8555.51, were strong enough to require explicit inclusion 
in the calculation. The scaling factor g(x) had to be applied 
to them, as well as A8542. Otherwise, with the Ca/N 

tot as- 
sumed for the deep photosphere, they would have been huge. 

Fig. [TU] shows the third member of the IRT, A8662. No 
adjustment of the parameters from those used to fit the 
A8542 line has been made. 
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Figure 10. The Call A8662-line, observed (solid gray) and cal- 
culated (solid black) with stratification. The calcium abundance 
and stratification parameters are identical to those used to fit the 
A8542-line. The model has log g = 4. 

7 STRONTIUM 

7.1 Weaker subordinate lines 

In addition to the resonance lines AA4077 and 4215, and 
the stronger subordinate lines AA4161 and 4305, we found 
five marginally unblended lines of Sr I and II, and determined 
provisional abundances from them. Results are shown in Ta- 
ble [T]for the 8600 K, log(g) — 3 model. The abundances are 
up to a factor of 3 lower if the log(g) — 4 model is used, but 
the agreement between the first and second spectra is still 
within 0.1 dex for the lower microturbulence. The weaker, 
Sri lines, which are less sensitive to microturbulence, mag- 
netic broadening, and surface gravity give mean abundances 
of —6.3 and —6.4 for the lower and higher values of £t re- 
spectively. 

We have not experimented extensively with stratifica- 
tion. The parameters used to get a fit for A4077 (see be- 
low) give unacceptably high abundances from the weak Sr I 
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Table 1. Abundances from weaker subordinate Sr lines: 8600 K, 
log(3) = 3. 



Spectrum 


A[A] 


W x [mA] 


log(Sr/N to t) 








£t = 1-3 ft = 3.5 


Sr II 


3380.69 


97 


-6.39 -7.70 


Sr II 


3474.89 


83 


-5.96 -7.23 


Sri 


6408.44 


23 


-6.21 -6.31 


Sri 


6504.00 


20 


-6.07 -6.15 


Sri 


6617.27 


18 


-5.87 -5.93 


Averages 






-6.10 -6.66 


log(Sr/Sr ) 






+3.02 +2.96 




4070 



4075 



4080 



4085 



Wavelength 



Figure 11. The Srn A4077-line, observed (solid gray) and cal- 
culated (solid black) with stratification. The parameters of the 
fit are: a = 1.7, logb = -6.70, d = -1, log(Sr/N to t) = -4.00. 
At log(r 50 oo) = 0.0, log(Sr/N to t) = -4.85. This is still consider- 
ably larger than the value obtained from the weak lines (Table [TJ: 
-6.10, with ft = 1.3. 



lines. The "standard" stratification parameters, a = 6.7, 
log b = —4.0, and d = +1, reduces the abundance difference 
between the first and second spectra though only by a little, 
and the overall abundances increase by more than a dex. 

We conclude that a reasonable mean abundance for the 
atmosphere of HR 710 is log(Sr/N to t) ~ —6.4. This is a fac- 
tor of 2.5 higher than the abundance found by Ryabchikova, 
et al. (1997) for the strontium star 7 Equ. It comports well 
with our assertion, based on Fig. [JJ that strontium is some- 
what stronger in HR 710 than 7 Equ. 



7.2 The resonance lines 

We were unable to fit the near wings of the stronger res- 
onance line, A4077, without moving the abundance jump 
deep into the photosphere. The adopted parameters given 
in Fig. 1111 lead to g(x) = 0.770 at the deepest layer of our 
atmosphere, log (T5000) = 1-4. The strontium fraction de- 
creases rapidly from the deepest layers, and levels to the 
value 6-(Sr/N tot ) = 2 ■ HT 11 by log (r 50 oo) = -5.4 

Fig. [12] shows the abundance stratification profiles of 
strontium and calcium. 

The weaker of the resonance lines, A4215, is shown in 
Fig. 1131 The calculation shows used the same fitting param- 
eters as those used for the stronger line. It is unclear just 




Figure 12. Stratified abundances of strontium (solid) and cal- 
cium (dashed) using parameters appropriate for the A4077 and 
K-line, respectively. Solar abundances are indicated. 




4205 



4210 



4215 
Wavelength 



4220 



4225 



Figure 13. The Srn A4215-line, observed (solid gray) and cal- 
culated (solid black) with stratification. The parameters are the 
same as those used in the fit for A4077. A considerably bet- 
ter fit (not shown) can be obtained from a slight modifica- 
tion of the parameters to: a = 1.7, log b = —6.70, d = — 1, 
log(Sr/N t ot) = -3.70. 



how meaningful this might be, since the blending is quite 
significant, and extends well into the core of the profile. 



7.3 The subordinate Srll lines 

Fig. ll4l shows a calculation of the subordinate line, A4161.80, 
5 P 2p i/2 ~ 6s 2 Si /2 . The line is in Srll Multiplet 3. Only the 
long-wavelength part of the profile is sufficiently unblended 
for use in fixing the stratification parameters. The calcula- 
tion shows the result of using the parameters from the best 
A4215-fit, which is poor. While we were able to improve the 
fit in the red wing somewhat, by trial-and-error adjusting 
the parameters of g(x), we never achieved a satisfactory fit. 
When the near wing, marked 'W in the figure, was fit, the 
calculated continuum (C) was below the observed one. 

A similar situation was found for the stronger line of 
Multiplet 3, A4305, though for this line it was the violet 
wing that proved most useful (cf. Fig. [TJ. 

In all likelihood, a satisfactory fit could be found with 
a more elaborate stratification model than the one fixed by 
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4160 



4165 



Wavelength 



Figure 14. The Srn A4162-line, observed (solid gray) and calcu- 
lated (solid black) with stratification. The parameters of the fit 
are: a = 1.7, log b = -7.0, d = -1, log(Sr/N to t) = -3.70, the 
same as for the A4077 resonance line. Other atomic lines have 
been omitted in the calculated profile for clarity. The near wing 
(W) and continuum (C) are noted. The observed violet wing con- 
tains a weak absorption, primarily due to Fe I (b). It is not strong 
enough to confuse the location of the near wing of the Sr II line. 



Table 2. Summary of stratification parameters 



A[A] 


log(El/N to t) 


a 


logb 


d 


Rem. 


Can 


3933 


-3.10 


6.7 


-7.00 


0.25 


Fig. 7 


3933 


-4.00 


6.7 


-6.00 


0.50 


logs = 4 


8498 


-2.12 


6.7 


-5.70 


0.25 


Fig. 8 


8542 


-1.52 


1.0 


-6.70 


-1.0 


Fig. 9 


8662 


-1.52 


1.0 


-6.70 


-1.0 


Fig. 10 


Sr II 


4077 


-4.00 


1.7 


-6.70 


-1.0 


Fig. 11 


4215 


-3.70 


1.2 


-6.70 


-1.0 


See text 


4162 


-3.70 


1.7 


-7.00 


-1.0 


Fig. 14 



the parameters of our p(a;)-function. We do not explore this 
here. 



8 STRATIFICATION AND HR 710 

8.1 Summary of stratification parameters 

Stratification parameters are collected in Table [2] The sec- 
ond entries for A3933 give a good fit to the core and wings 
(to ca. ±40 A) for a surface gravity \ogg = 4.0. There is 
no corresponding figure. The entries for Srn A4215 are the 
values mentioned in the caption to Fig. 1131 which give an 
improved fit. Note that (El/N to t ) refers to the deepest pho- 
tospheric layers. 



8.2 General remarks 

The near wings and cores of strong lines in magnetic CP 
stars cannot be fit with traditional models. Non-LTE is 
surely relevant, but it is most unlikely to account for all 
of the discrepancies between theory and observation. Mod- 
ern workers have returned to the notion of "anomalous at- 
mospheres, or AA," a term coined by Underhill (1966). 
Kochukhov, Bagnulo, & Barklem (2002) showed that the 
Balmer-line profiles could be fit with an ad hoc modifi- 
cation of the temperature distribution in the high photo- 
sphere. Ryabchikova and colleagues (cf. Ryabchikova 2005) 
have used stratified abundance models to explain both the 
cores and near- wings of strong lines of Ca n and anomalous 
excitation and ionization, often seen in magnetic CP stars. 

One might expect that the same stratification model 
would fit all of the Can lines, H and K, and the IRT. This 
was not the case for HR 710. We must note that our strat- 
ification model is somewhat more constrained, by Equ. 1, 
than those used by Ryabchikova and coworkers. 

The cores of the resonance lines of Srn, AA4077 and 
4215 exhibit similar peculiarities to the K-line. Vastly im- 
proved fits can be obtained using stratified models, though 
slightly different stratifications are required. The subordi- 
nate lines, AA4162 and 4305 require stratification or other 
modifications to fit the wings. The stratification parameters 
required to fit A4077, however, make very little difference in 
the wings of the two subordinate lines. 

Ionization and excitation anomalies are present in the 
HR 710 spectrum. The Til and II ionization and excitation 
imbalance deserves further study. The rich Fe I and II spectra 
were in satisfactory agreement without the need for strati- 
fication, provided we employed the significantly low-gravity 
model with \ogg = 3.0. 

Generally, the line spectrum in the red and infrared was 
fit better by a higher-gravity model, and the relevant illus- 
trations shown are based on that model. We cannot rule 
out the relevance of abundance patches for our difficulties 
in fitting all of the Ca II or the Sr II lines with a single strat- 
ification model. 



9 HR 710AA AND ITS COMPANION HR 
710AB 

The mass function for an SB1 system is f(m) = 
m 3 sin 3 (i)/(m + M) 2 Mq, where M and m are the masses 
of the primary and secondary, respectively. This /(m) was 
well determined by Bonsack, and has the value 1.777 ■ 10~ 3 . 

For a given m and f(m), the orbital inclinations de- 
crease as M —> m. Kitamura (1980) noted that the inclina- 
tions of short-period Ap SB1 stars tend to be small. Such 
near pole-on views of SB1 systems would require a sublumi- 
nous star, since the companion spectrum is unseen. Bonsack 
suggested a white dwarf for HR 710. This model explains the 
failure of the magnetic field to show an obvious reversal. 

If we assume the 3-day orbital period is also the period 
of rotation of the primary, the observed v ■ sini depends 
directly on the radius of the star and the inclination of its 
rotational axis. 

The equatorial velocity of a star is 

50.56 • R 

v cq = 5 , (2) 
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Table 3. Secondary mass m for various primary masses M and 
v ■ sin i. The preferred values are in bold face. 



M/Mq 


R/R Q 


i (deg) 




V ■ sini 










8 


10 


12 


2.00 


1.74 


29 


0.89 


0.68 


0.56 


1.84 


1.63 


28 


0.79 


0.62 


0.48 


1.60 


1.46 


25 


0.64 


0.49 


0.40 



where 7? is in solar radii and P in days. 

Values of v-sin i have been published by Bonsack (1981) , 
and Preston (1971): 7.7, and 13 kms -1 respectively. We have 
made a new determination from our UVES spectra, using 
the magnetic null lines AA3850, 4006, 4065, and 7090, all of 
Fe I. 

The observations were (Fourier) unfiltered, but they 
were interpolated to every 0.02 A, for compatibility with 
Michigan software. We assumed a Gaussian 1/e- width of 
0.04 A for the instrumental profile in the blue-violet. This 
allows for a small degredation of the resolution; a resolving 
power, 1Z of 80 000 would give ca. 0.03 A for the 1/e- width, 
where we use A A = X/TZ as the FWHM of the instrumen- 
tal profile. For A7090, using TZ = 110 000, the 1/e width 
of the pure instrumental profile is 0.045 A, and we assumed 
0.055 A. 

Synthetic spectra were fitted to regions containing these 
lines spectra, and the best fits determined by trial and error. 
We find that a value of 10 kms" 1 gives the best overall fit. 
This is near the mean of the Preston and Bonsack's pub- 
lished values. However because of blending, and subjective 
judgements about a fit, the uncertainty of v ■ sini is still 
between 1 and 2 kms" 1 . 

We prefer M/M© = 1.84 and R/R Q = 1.63 from tables 
in Drilling and Landolt (1999), for a main sequence A-star 
with an effective temperature of 8600 K. The values agree 
well with Hubrig, North, and Scholler (2007). 

It is now straightforward to get the inclination from 
Eq. (2), if we assume a 3-day period. From the inclination, 
the secondary mass, m, will follow from the mass function, 
given a primary mass, M. In order to provide some indica- 
tion of the overall uncertainties, we give the secondary mass 
in Table [3] for three assumed masses (and corresponding 
radaii), and three values of v ■ sini. We assume the relation 
between mass and radius for main sequence stars given by 
Drilling and Landolt. 

Our preferred values give m — 0.62 Mq, conveniently 
near the most probable observational value for white dwarf 
masses according to Weidemann (1990). This supports the 
assertion of Czesla and Schmitt (2007) that the secondary 
of HR 710Aa "is clearly a very plausible candidate for the 
observed X-ray emission...". These authors discuss a mag- 
netically confined wind-shock model. 

Bonsack (1981) also discussed a 15.88-day magnetic pe- 
riod. With R — 1.63 Rq, this period leads to an equatorial 
velocity of 5.2 kms -1 , incompatible with the adopted v-sini. 
To reconcile this period with the observed v ■ sin i, we would 
need to assume a stellar radius nearly double that used. This 
is not compatible with the luminosity of the star. 



10 REMARKS AND CAVEATS 

The present work on HR 710 shows that it is difficult to 
explain the profiles of strong Ca II and Sr II lines with a sim- 
ple abundance-stratified model. Overall, we strongly support 
the stratification models, but find indications that refine- 
ments are necessary. 

Different stratification parameters were required for the 
K-line and the IRT of Can, while Srll lines also required 
adjustments of the parameters for individual lines. 

It remains to be seen if theoretical work can explain 
the large elemental abundances of the best-fit stratification 
models. These require very high abundances in the deep 
photosphere in order to explain the observed, broad wings. 
In the case of calcium, traditional diffusion calculations have 
explained low photospheric abundances (Am stars). 

Among a number of uncertainties in the present dis- 
cussion, we note that our stratification model (g(x)) may 
(a) not have explored the optimum domains of parameter 
space, or (b) have been too simplistic. Regarding the latter 
point, we note the sophisticated approach of Kochukhov, et 
al. (2006). 

Eventually, it will prove necessary to approach the mod- 
eling of the spectra of these stars with improved meth- 
ods. Useful modifications of our procedures might include 
changing the T — r relation in the high photosphere. This 
has already been shown to give fits to the Balmer lines 
(Kochukhov, Bagnulo, & Barklem 2002), and can give nibs 
on the K-line (Cowley, Hubrig, & Kamp 2006). In the last 
reference, it was also pointed out that sharp K-line nibs 
could result if the gas pressure in the upper atmosphere 
were lowered. 

Calculations done in LTE must be considered provi- 
sional until supported by work done with non-LTE. This is 
especially true for the upper atmospheres, where the total 
radiative flux is very nearly constant, independent of local 
conditions. The core of the Ca K-line is formed highest of 
lines considered in this work, and Cowley, Hubrig & Kamp 
showed that LTE calculations were only slightly modified by 
non-LTE. 

The greatest challenge of the observations of strong-line 
profiles of CP stars is to the theoretician. It is straightfor- 
ward to make ad hoc adjustments to a model until a profile 
is fit. It is quite another to justify such modifications with 
rigorous calculations. 
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APPENDIX A: CALCULATION OF HIGHER 
PASCHEN LINES 

Calculations near series limits need to account for overlap- 
ping lines, as well as the "dissolution" of some discrete up- 
per levels by various perturbations from the plasma bath 
of the emitting and absorbing atoms. We have implemented 
the scheme outlined in Cowley (2000) in conjunction with a 
probability function W(n) for the existence of discrete upper 
levels. The specific routine was from Barklem (2004). The 
basic assumption is that we can calculate the total opacity 
due to both bound-bound and quasi-continuous absorption 
by considering only a wavelength region near a given bound- 
bound transition. These regions extend from the halfway 
points between transitions to level n — 1 and n, and from n 
to n+1. We shall refer to these regions as the n-local regions. 
Following the reference cited, we define an n-dependent cross 
section, such that for a Paschen line 

/ CT 3 nrfAA = ^-^\ 2 n fzn- 

J me* 

The integral here is assumed to extend only over the n-local 
region. The cross sections so defined are in excellent agree- 
ment with those found in TOPbase (Cunto, et al. 1993) for 
photon energies below the ionization limit. 

Line absorption in the n-local region is calculated in the 
usual way, but with the product Nsf^n diminished by the oc- 
cupation probability W(n). However, the line profile is trun- 
cated at each end of the n-local region. The ordinary con- 
tinuous opacity is then augmented by an amount necessary 
to bring the total absorption into agreement with the above 
equation. In this way the total absorption merges smoothly 
with the known photoelectron absorption cross section at 
the series limit, as has been known for many years (Sugiura 
1927). 



